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ABSTRACT 

Recent results indicate that the grain size and crystallinity inferred from observations of silicate 
features may be correlated with spectral type of the central star and/or disk geometry. In this paper, 
we show that grain size, as probed by the 10 /xm silicate feature peak-to-continuum and 11.3-to- 
9.8 /xm flux ratios, is inversely proportional to log L*. These trends can be understood using a simple 
two-layer disk model for passive irradiated flaring disks, CGPLUS. We find that the radius, i?io, of 
the 10 /ira silicate emission zone in the disk goes as (L^/Lq) ' 56 , with slight variations depending on 
disk geometry (flaring angle, inner disk radius). The observed correlations, combined with simulated 
emission spectra of olivine and pyroxene mixtures, imply a grain size dependence on luminosity. 
Combined with the fact that Rio is smaller for less luminous stars, this implies that the apparent 
grain size of the emitting dust is larger for low-luminosity sources. In contrast, our models suggest 
that the crystallinity is only marginally affected, because for increasing luminosity, the zone for thermal 
annealing (assumed to be at T > 800 K) is enlarged by roughly the same factor as the silicate emission 
zone. The observed crystallinity is affected by disk geometry, however, with increased crystallinity 
in flat disks. The apparent crystallinity may also increase with grain growth due to a corresponding 
increase in contrast between crystalline and amorphous silicate emission bands. 

Subject headings: circumstellar matter — stars:pre-main sequence-infrared: protoplanetary disks: lines 
and bands — stars: formation — solar system: formation 



1. INTRODUCTION 

Silicate emission features from circumstellar disks have 
now been observed toward large samples of young stars, 
ranging from sub-stellar mass brown dwarfs to Her- 
big Ae/Be (hereafter HAEBE) stars up to 10 times 
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the mass of the sun (see, e -g.. iBouwman et alJ 
van Boekel et alJl2003l l2005t IKessler-Silacci et all 
2006t lApai et all 120051 ) . The most commonly observed 
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silicates are amorphous olivines and pyroxenes and 
their crystalline, Mg-rich forms, enstatite (MgSiOa) and 
forsterite (Mg2Si04). For amorphous olivine and pyrox- 
ene, the features due to distinct stretching or bending 
modes of the silicon-oxygen bonds are blended into two 
broad features near 9.8 ^m (stretch) and 18 ftm (bend). 
In addition to providing information about the grain 
composition and crystallinity, sili cate features have been 
used as a probe of the grai n size (|van Boekel et al.ll2003t 
IKessler-Silacci et"aTI l2006t ). Since the silicate emission 
features arise from the optically thin surface layer of the 
disks, and larger grains are expected to settle to the disk 
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midplane, the grain size and crystallinity indicated by 
these features applies primarily to the smallest grains in 
the disk. 

Studies of silic ate emission fe atures around HAEBE 
stars and comets (jWood cn 2002), indicated that mature 
disks and comets had silicate emission features character- 
istic of up to 30% crystalline silicates, while much lower 
crystallinity was obse rved toward many young disks 
(IBouwman et al.ll2001h. embedd ed objects (jBowev et al.1 
Il998t IWhittet fc Tielensl Il997fl and through the ISM 
toward the galactic center (|Kemper et al.1 12004. 2005). 
This suggests an evolutionary trend in which amorphous 
silicates from the ISM are crystallized within circumstel- 
lar disks. Indeed, the types of silicate features observed 
seem to support this scenario. Van Boekel et al. (2005) 
found that sources showing large mass fractions of crys- 
talline silicates also possessed higher mass fractions of 
large grains, suggesting that crystallinity is somehow re- 
lated to grain coagulation in disks - this is not necessarily 
a causal relationship. For disks around sub-stellar mass 
stars, the spectra of flat disks, in which dust coagulation 
and settling have resulted in smaller flaring angles, of- 
ten show higher crystallinit y and larger grai n sizes than 
thos e of more flared disk s (lApai et al.ll2005l ). Addition- 
ally, IKessler-Silacci et al.1 (|2006f) show that the grain size 
indicated by the 10 fim feature is larger for disks around 
M-stars than for disks around A/B-stars, which they in- 
terpret as an indication that the silicate feature probes 
different regions of the disk depending on stellar lumi- 
nosity. 

In order to examine the physical processes responsible 
for these trends, we use a 2-layer disk model to trace the 
silicate emission zones for a range of stellar masses from 
0.003-7 M . The relationship between observed silicate 
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the flattening of the ob- 
by the 11.3-to-9.8 /jm flux 



emission feature characteristics and stellar luminosity is 
described in §2. The disk models are described in §3. 
Finally, in §4, the results and implications are discussed. 

2. SILICATE EMISSION FEATURES, GRAIN SIZE AND 
STELLAR LUMINOSITY 

The connection betwee n silicate emission s treng th and 
shape was first noted by Ivan Boekel et al] ((2003) for a 
sample of HAEBE stars, and has now been seen for 
a number of sources, including HAEBE stars, T Tauri 
stars and, more recently, a small sample of brown dwarfs 
(IPrzvgodda et aUl2003t IKessler-Silacci et al.ll2005l . 120061 : 
lApai et all 120051 ) . These results are shown in Figure [IJ 
with the addition of ~40 new sources more recently ob- 
served by Olofsson et al. (in prep.) as part of the 
"From Molecular Cores to Planet-Forming D isks (c2d)" 
Spitzer Legacy program (| Evans et al.l |2003[). including 
9 brown dwarf candidate s selected from lAllers et al.1 
(I2006h. As was sho wnhi lvan Boekel et all (|2003f ) and 
IKessler-Silacd et al.1 (|2005[ ) 
served shape, indicated 
(S11.3/ Sg.s), and decrease in the observed strength 
i^peaD °f the silicate features can be tied to increases 
in the amorphous olivine/pyroxene grain sizes, such that 
grain size increases from the bottom right to the top 
left of Figure [T] Crystalline forsterite and/or PAH 
emission can result in increased flux near 11.3 /zm 
and therefore cause devia tion from this trend (see also 
iKessler-Silacci et al] I2006T) . In particular, spectra with 

(Sii.s/Sg.sV^r > 1-0-1.1 (approximately l/5 th - 
1 /3 rd of the sample shown in Figure [T|) cannot be ex- 
plained with models of pure amorphous silicates. Over- 
laid in Figure [1] are parameters calculated from modeled 
spectra of pure olivine and olivine-pyroxene (75%:25%) 
mixtures for filled homogeneous spheres (usin g Mie the- 
ory) or distribution of hollow spheres (DHS (Mi n et all 
|2005[ ) , representing porous grains) . The olivine-pyroxene 
mixtures correspond well with the observed trend. 

Figure [T] shows a trend in which the disks around the 
lowest luminosity brown dwarfs have larger grain sizes 
than those around T Tauri or HAEBE stars. This is 
made more apparent in Figure [2^,, where we plot the 
grain size parameter versus the logarithm of the stel- 
lar luminosity for the same sample as shown in Fig- 
ure [TJ Indeed, although the exact relationship between 
grain size and (Sn.3/ 'Sg.s) / 'S pe £™ may vary depending 
on the grain composition/shape (Figure O)), there ap- 
pears to be an inverse correlation between the grain size 
parameter and the stellar luminosity. Brown dwarfs have 
larger (Sj.1.3 / Sg.s) / S^J^™ than T Tauri stars, which have 
larger (Sii.a/Sg^/S^™ than HAEBE stars, and al- 
though there is a large degree of scatter, the correlation 
is significant, with a >97% probability that x and y are 
correlated and all data lyin g within 2a of the fitted trend. 
IKessler-Silacci et al.l (|2006h postulated that the observed 
relationship between grain size and spectral type could 
be explained if the silicate features probed different dis- 
tances from the star, probing smaller grains located far- 
ther out in the disks around more luminous stars. 

Using the relation between (Sn.3/S9.g)/S pe £™ and 
log L+ from observations in Figure [2^ and the 
{S11.3/ Sg.s) /^peafe" an0 - g ram s i ze from modeled spectra 
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FlG. 1. — Silicate 10 fj,vn feature strength vs. 
squares are brown dwarfs with M* < 75 Mj up . 
icate features show up toward the right of the plot and fiat 
features show up toward the left. Blue squares are low-mass 
stars with M+ > 75 Mj up . Red triangles are HAEBE stars. 
Filled symbols are from the S p itzer c 2d sample. Open sym - 
bols are from Ivan Boe kel et all 12001); I Przvgodda ~ al. (2003); 
IKessler-Silacci et al.l (|2005l ); lApai et al J 120051) ■ The black line 
denotes a fit to the points with y < 1.1, y = — (0.18±0.02)x + 
(1.23±0.03), with r = -0.77, indicating a probability of >99.05% 
that x and y are correlated. Models of pyroxene and olivine mix- 
tures, treating the dust as filled homogeneous spheres (Mie) or dis- 
tribu tions of hollow spheres (DHS ) with grain sizes from 0.1-20 fim, 
from Kessler-Silacci ct al. (2006) are overlaid as dashed/dotted 
lines. 

in Figure ^jp, one can determine a relationship between 
grain size and luminosity, shown in Figure [2J;. The fea- 
ture strength and shape for a particular grain size can be 
quite differ ent for different silicate co mpositions /shapes 
(Figure l2"b lKessler-Silacci et al.|[2006l ). This in turn af- 
fects the precise relationship between grain size and stel- 
lar luminosity (Figure^). In particular, the DHS models 
show stronger features than the Mie models for the same 
grain size, and thus for low luminosity stars, the DHS 
models require grain sizes that are much larger (~20 /im) 
to produce the weak features. Nevertheless, the conclu- 
sion that apparent grain size is inversely proportional to 
luminosity is robust for all reasonable grain models. 

3. MODELS OF THE SILICATE EMISSION ZONE 

In order to understand possible physical explanations 
for the correlation shown in Figured, we model the disks 
around a set of stars ranging from -10" 4 -10 4 L Q to bet- 
ter identify how the silicate emission zone — the portion 
of the disk contributing to the 10 /mi feature — changes 
as a function of the stellar lum i nosity . We use the 2- 
layer model of iDullemond et al] (l200ll) (CGP LUS). which 
is based on work by iChiang fc Goldreichl (l997) and 
IChiang et all (|2001h . This model assumes vertical hy- 
drostatic equilibrium, resulting in a flared disk structure. 
Directly illuminated dust grains in the optically thin sur- 
face layer re-radiate half of the energy absorbed from the 
central star down into the disk interior and thus regulate 
the interior disk temperature. The other half is emitted 
away from the disk and can be observed as optically thin 
emission from the disk (e.g., emission features). We em- 
ploy a new version of CGPLUS in which the originally 
vertical blackbody puffed-up inner dust wall model is re- 
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Fig. 2. — (a) Grain size parameter vs. stellar luminosity. Symbols are as in Figure [T] The black line denotes the correlation, y = 
— (0.05±0.02)a; + (0.61±0.02), with r = —0.23, indicating a probability of >97% that x and y are correlated. (6) Models of grain mixtures 
(as in Figure [TJ relating the grain size to {Sxi.3/Sg,a)/S^^. (c) Models relating log(L*/LQ) to grain size, using the relation derived 
from (a) applied to the models in (b). Line types are the same as in (6). 



placed by a smoothly rising disk height from H = at the 
dust evaporation radius to H = U kT m idR 3 / fimpGMi, at 
a ~20% larger radius, where T m ^ is the midplane tem- 
perature at radius R, \i is the gas mea n molecular weight, 
and m p is the mass of a proton (see ICalvet et al.l |2002| 
for a description of this method to mimic a rounded-off 
inner dust rim) . The location of the inner rim of the disk 
is determined by the dust sublimation temperature (typ- 
ically 900-1 600 K for silicates, a ssumed here to be 1500 K 
as found bv lDuschl et al.lfl996l ). The inner radius is set 
to R+ if the dust sublimation tem perature is not reached . 
Opacities used in CGPLUS are from lLaor fe Draind (l993), 
for mixtures of 0.1 /an spherical amorphous olivine and 
crystalline graphite grains. 

To enable comparison between the models, the disk 
mass is assumed to be a simple function of the stel- 
lar mass, M disk = 0.03 M« . Using M disk (xM? 
might be a better ass umption (jHueso fe GuillotJ 120051 : 
iDullemond et al.1l2006[ ). but has a very small effect (an 
offset of < 0.02 AU) on the radius of the silicate emis- 
sion zone. Stellar luminosities, masses and surface gravi- 
ties are re lated using th e pre-main-se quence evolutionary 
tracks of iBaraffe etaLl (11998L 12003(1 for = 0.0002- 
0.5 L Q and lSiess et all 1)20001 ) for L+ = 0.1-2000 L . To 
check consistency, both models are used from = 0.1- 
0.5 Lq, and we find no discontinuities in the overlap re- 
gion (see Figured]). The outer radius of the disk cor- 
responds to where gas would no longer be bound to 
the star (0.2r o , where r g = GM*/j,m p /kT; see, e.g., 
IDullemond et al.l l2~006a). The disk surface density is de- 
scribed as £ = E (i?/1 AU)' 9 , where R is the radial dis- 
tance from the star, So is a constant set by the disk mass, 
and j3 = — 1 for all disks. Due to the fact that we ex- 
tend the disk out to the radius where the gas becomes 
unbound, the disk in these outer regions will be very 
geometrically thick. This makes it difficult to explore 
inclination effects using the simple two-layer model, and 
true 2-D/3-D disk models are beyond the scope of this 
paper. Instead we always assume that the disk is face-on 
(i = 0°) in our models. However, variations in the in- 
clination angle could contribute to the scatter shown in 
Figure [2^, with higher inclinations resulting in decreased 



10 feature strength. 

CGPLUS produces both spectral energy distributions 
and images at each wavelength for each model. In order 
to evaluate the 10 firm emission zone, we calculate the 
radially integrated cumulative flux of the 10 fiisi emis- 
sion image. The 10 fim emission zone is then defined as 
the radial zone where the cumulative flux reaches 20%- 
80% of the total 10 /xm emission. Figure EK, shows the 
10 /jm emission zone (yellow solid region) for the result- 
ing models. The location of the zone increases roughly 
as Rio = 0.35 AU (L*/L Q ) - 56 . This relationship can 
be easily approximated under the assumption of equi- 
librium between the absor bed and emitted energy (e.g., 
iBackman fc Parescdll993l equations 1 and 2), with the 
simplifying assumptions that 1) the emission is coming 
from an optically thin region of the disk, 2) only a nar- 
row range of temperatures is probed at 10 /im, and 3) 
that the grains absorb to first order as blackbodies. The 
resulting distance between the illuminating star and the 
emitting grains is R oc L 5 , which is very similar to the 
relation found here. 

The silicate emission therefore arises from very dif- 
ferent regions in disks around brown dwarfs (-Rio < 
0.001 - 0.1 AU) than in disks around HAEBE stars 
(-Rio > 0.5 — 50 AU). It is likely that the conditions 
in these disks, and thus properties of the silicates, are 
very different at these radii. We will now use some sim- 
ple approximations to qualitatively indicate how radial 
conditions can affect the silicate properties (grain size, 
crystallinity) that are probed with the 10 /im features. 

3.1. Effects of Radial Conditions on Grain Size 

Figure shows that the grain sizes interpreted from 
observations of the 10 /jm feature appear to be related 
to th e stellar luminosity (see also rKcsslcr-S ilacci et al.l 
2006). To simplify, we will attempt to explain the ob- 
served trends using Ansatz that the grain size distribu- 
tion with radius is the same in all disks, and that only the 
silicate emission region changes. We are aware that this 
may not be quite realistic, but co agulation is at present 
insufficiently understood (see, e.g.. lDominik et al. 1120061 
for review) and we therefore wish to study - in isolation - 
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Fig. 3. — Radius probed by 10 fim feature as a function of stellar luminosity, (a) The yellow solid region shows the radii contributing 
20—80% of the 10 fim emission for flared disks, (b) The additional blue horizontal striped and green vertical striped regions show the same 
for flat disks and flared disks with inner radii 3 times larger than the dust sublimation radius, respectively. In both panels, the black dotted 
and red dashed lines show the radii corresponding to T = 1500 K and T = 800 K, above which temperature silicates will sublimate and 
crystallize via annealing, respectively. The vertical black dashed lines show the luminosities corresponding to grain sizes of 3, 4, 6, and 
10 fim for the grain size vs. stellar luminosity trend in Figure [2J:. The horizontal black dashed lines show the radii corresponding to those 
luminosities, using the relationship shown by the filled regions, log(i?) = —0.45 + 0.561og(L*/LQ). 



the effect of the effect of the radial location of the silicate 
emission zone on the observed grain sizes. 

The grain size distribution vs. luminosity deduced from 
observations and shown in Figure^, has been used to in- 
fer the effect of radial conditions on grain size. The rela- 
tion for DHS olivene-pyroxene mixture has been adopted 
since this model best represents the trend shown in Fig- 
ured Using the modeled variation, the grain sizes of 3, 4, 
6 and 10 /im would correspond to the luminosities shown 
as dashed vertical lines in Figure [3] Using the modeled 
variation of the location of the silicate emission zone with 
stellar luminosity shown by the filled regions, we then 
calculate the radii corresponding to those luminosities 
(shown as dashed horizontal lines). This results in grain 
sizes that decrease as a function of disk radius, which 
seems reasonable since coagulation and dust settling pro- 
cesses are believed to occur more quick ly at smaller radii 
(see, e.g. iDullemond fc Dominikl I2005D . Thus this "re- 
verse modeling" shows that the observation that later 
type stars appear to have larger grains (as deduced from 
their spectra) might be due, at least in part, to the fact 
that grains are expected to be bigger at smaller radial 
position in the disk. 

The example presented above is indicative of how stel- 
lar luminosity, and the corresponding location of the 10 
/im silicate emission zone, might affect the grain size de- 
duced from silicate emission features. However simple, 
this example illustrates that the location of the silicate 
emission zone is dependent on stellar luminosity. This 
effect cannot be ignored when one seeks to interpret 
the grain sizes derived from silicate emission features in 
terms of grain growth in the disk as a whole. 

3.2. Effects of Radial Conditions on Crystallinity 

Although crystalline silicates have been observed in 
disks around HAEBE stars, T Tauri stars, and brown 
dwarfs, the mechanism for conversion from the primarily 
amorphous ISM silicates is not completely understood. 
Spatially resol ved mid-IR interferomet ry of disks around 
HAEBE stars (jvan Boekel et all 120041 ) indicate that the 
crystallinity is much higher at small radii (< 2 AU) than 



in the outer disk. This result supports the theory of 
crystallization via annealing at high temperatures in the 
inner disk, followed by transport to larger radii. In this 
case, the amount of crystalline silicates present in the 
10 /im emission zone is a function of both the radial lo- 
cation of the zone and the timescale for radial mixing. 

To understand the effects of probing different radii on 
the observed 10 /im crystallinity, we plot in Figure [3ji 
the radii corresponding to T s = 800 K, the minimum 
temperature required for ann ealing amorphous olivine 
(|Duschl et all 119961 lGaillll998h . If crystallization is oc- 
curing primarily through annealing, then we would ex- 
pect dust to be crystallized at radii equal to or smaller 
than the dashed red line in Figure^. As the 10 /im emis- 
sion arises from the disk surface layer, we use only the 
temperature of the surface layer to calculate the region. 
For simplicity, the default version of the CGPLUS model, 
with a single small grain size throughout the disk, is used 
in these calculations. 

Although the 10 /im emission zone moves to larger radii 
for HAEBE stars with respect to brown dwarfs, the crys- 
tallization zone also expands to include larger radii; more 
luminous stars can heat dust farther out in the disk. 
Thus, we expect the percentage of crystalline silicates 
emitting at 10 /im to be nearly independent of the lu- 
minosity of the central object. There does appear to be 
a small increase in R cr yst / Rio with increasing stellar lu- 
minosity, however, such that we might expect to observe 
more crystalline silicates from HAEBE star disks than 
brown dwarf disks. In contrast, observations show that 
crystalli nities can be qu ite large in disks around brown 
dwarfs (jApai et al.ll2005t Merin et al. 2006, submitted). 
This difference could be explained by accretion, which 
would result in viscous heating of the midplane of the 
inner disk, thus leading to additional crystallization at 
smaller radii. Additionally, the degree of crystallinity in 
these disks could be affected by the presence of additional 
crystallization methods or the effectiveness of radial mix- 
ing, which are not considered here. Overall, the model 
predicts nearly constant crystallinity as a function of lu- 
minosity. 
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3.3. Effects of Disk Geometry 

Because disk geometry affects the incidence of stellar 
light onto the surface of the disk, it should also affect 
the absorption and resulting emission of this light via the 
10 /xm feature. The effects of disk geometry are shown 
in Figure [3)d. The yellow solid region shows the 10 fim 
emission zone for models of flared disks, as in Figure^. 

The blue horizontal striped region shows the 10 /mi 
emission zone as calculated for models of flat disks, possi- 
bly resulting from dust settling, with otherwise the same 
parameters as the flared disks. The flaring angle of the 
disk determines the angle at which the stellar radiation 
penetrates the disk surface. Since the radiation only 
penetrates to t — 1 (along the line of sight), the true 
(vertical) depth decreases with the incidence angle cf> as 
sin(</>). So for a flat disk (with a steep drop-off of (j) 
with radius), the depth at which the stellar radiation 
heats the disk decreases very steeply with radius. This 
means that the amount of matter contributing to the 
10 /xm feature also decreases much faster as a function 
of radius for a flat disk than for a flared disk. So the 
10 /xm feature from a flat disk is dominated by emission 
from shorter radii, while the feature from a flared disk 
includes also emission that is spread out over a larger 
portion of the disk. With our prescriptions for grain size 
and crystallinity, this means that the 10 /xm emission 
from flat disks would probe predominantly smaller radii, 
wi th larger and slig htly more crystalline silicates, as seen 
in lApai et all (|2005l ). 

The green vertical striped region in Figure [3Jd shows 
the 10 /im emission zone for models of flared disks with 
inner radii that are larger than the dust sublimation ra- 
dius. By moving the inner radius of the disk outward, we 
simulate the effect of clearing the inner disk, e.g., by a 
planet or photoevaporation. Models of this type success- 
fully reproduce spectral energy distributions observed for 
several T Tauri s tar disks, such as GM Aur, DM Tau 
and CoKu Tau 4 (ICalvet et al.ll2005t iForrest et aLll2004t 
iD'Alessio et al.ll2005l) . As" seen in Figure [3)3, the addition 
of an inner hole in these models means that the 10 /im 
emission is now probing larger radii (especially for low 
luminosity stars). 

As shown in Figure [3Jd, the disk geometry (flaring an- 
gle, inner holes) can affect the radial location of the 
10 /xm emission zone. The presence of a variety of disk 
geometries could contribute to the observed scatter in 
the trend noted in Figure [21 while preserving the overall 
correlation with L*, as observed. 

4. DISCUSSION AND IMPLICATIONS 

In this paper we have demonstrated that the region 
of a protoplanetary disk that is probed by 10 /xm sil- 
icate emission features is related to the luminosity of 
the central object; specifically, the 10 /xm emission zone 
is at larger radii for disks around more luminous stars, 
Rio = 0.35 AU (L*/L ) - 56 . As a result, the grain sizes 
inferred from observations of silicate emission features 
may reflect the location of the silicate emission zone more 
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than the properties of the bulk silicate dust in the ob- 
served disk. If the grain size is a function of the ra- 
dial position in the disk (as shown here by extrapolating 
from the observed feature shape/strength vs. luminosity 
and the modeled grain size vs. shape/strength relations), 
then the corr elation of grain size w i th ste llar type seen in 
the sample of lKessler-Silacci et al.1 (|2006| ) and the gener- 
ally large grain sizes inferred fr om observat i ons o f dust 
around very low mass stars by lApai et all (|2005T ) may 
be related directly to the location of the 10 /xm emission 
zone in these disks. Additionally, the 10 /jm emission 
zone can vary with disk geometry (flaring angle, inner 
disk radius), which may result in deviations from the 
trend for individual objects. 

The amount of crystalline silicates contributing to the 
observed 10 /xm feature, if dependent primarily on the 
disk temperature, would not be strongly affected by stel- 
lar luminosity, as both the crystallization region and the 
10 /xm emission zone increase in size with increasing L+. 
However, if the grain size is a function of radial position 
in the disk, then the apparent crystallinity inferred from 
the 10 /xm feature may be affected, due to the fact that 
the contribution to the spectrum from amorphous sili- 
cates diminishes more strongly with grain growth (0.1- 
10 /xm) than does that f rom crystalline silicates (c.f., 
iKessler-Silacci et al.ll2006l Figures 6-8). This effect may 
contrib ute to the large cr ystallinities seen toward brown 
dwarfs (|Apai et all 120051 : Merin et al. 2006, submitted) 
and the restriction of significant crystalline mass frac- 
tions (>15%) to those HAEBE disks that also possess 
>80% big (1.5 /xm) grains by mass (jvan Boekel et al.l 
[20051) . 

The relationship of the location of the 10 /xm emission 
zone with stellar luminosity and disk geometry demon- 
strated here will also affect studies of correlations of 
silicate properties (as deduced from silicate emission 
features) with stellar parameters, age, and disk evolu- 
tion. In particular, relationships of grain properties with 
disk/stellar evolution or age are difficult as the stellar lu- 
minosity and disk geometry are expected to change with 
evolutionary state of the system and both are shown here 
to affect the location of the 10 /xm emission zone. This 
study indicates the necessity for careful modeling of disks 
in order to deconvolve the effects of variations in the 
10 /xm emission zone and true changes in silicate dust 
properties. 
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